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Abstract: Previously it was shown that there exists a class of viscous cosmological models which 
violate the dominant energy condition for a limited amount of time after which they are smoothly 
connected to the ordinary radiation era (which preserves the dominant energy conditions). This 
violation of the dominant energy condition at an early cosmological epoch may influence the slopes 
of energy spectra of relic gravitons that might be of experimental relevance. However, the bulk 
viscosity coefficient of these cosmologies became negative during the ordinary radiation era, and 
then the entropy of the sources driving the geometry decreases with time. 

We show that in the presence of viscous sources with a linear barotropic equation of state p — yp 
we get viscous cosmological models with positive bulk viscous stress during all their evolution, and 
hence the matter entropy increases with the expansion time. In other words, in the framework 
of viscous cosmologies, there exist isotropic models compatible with the standard second law of 
thermodynamics which also may inffuence the slopes of energy spectra of relic gravitons. 

PACS numbers: 98.80.Cq, 04.30.Nk, 98.70.Vc 



Our universe can be viewed as containing a sea of 
stochastically distributed gravitational waves of primor- 
dial origin. Among all observational cosmological ev- 
idences (present and future), primordial gravitational 
waves should have a sufficiently enlightened character in 
order to better understand the very early universe. The 
gravitational waves of cosmological origin are nothing but 
squeezed states of many gravitons produced from the vac- 
uum fluctuations of the background metric. A qualita- 
tive analysis can be performed in the context of different 
physical frameworks, since all models for the very early 
universe predict the formation of stochastic gravitational 
wave backgrounds. As examples we can mention infla- 
tionary quintessential models inflationary models in 
Brans-Dicke theory of gravity cosmological models 
in the Brane- world scenario Q, and superstring theo- 
ries . The shape of the stochastic graviton background 
spectrum is affected by the variations of the background 
dynamics. 

In this context, Giovannini j5j has considered the 
interesting possibility of constructing flat Friedmann- 
Robertson- Walker (FRW) cosmologies endowed with a 
bulk viscous stress which induces a violation of the dom- 
inant energy condition (DEC) for a limited amount of 
time at an early cosmological epoch. This kind of cos- 



mological models may be connected to some of the re- 
cent remarks of Grishchuk |^ concerning the detectabil- 
ity of stochastic gravitational wave background by forth- 
coming interferometric detectors, such as LIGO, VIRGO, 
GEO600, LISA 0. Effectively, bulk viscous dissipative 
processes may influence the slopes of the energy spectra 
of relic gravitons (generated at the time of violation of the 
DEC) producing an increasing with frequency in a calcu- 
lable way. These slopes are crucially related to the sign 
of the p + p, where p and p are, respectively, the energy 
density and the pressure density of the cosmic fluid. The 
requirement that one wants expanding and inflationary 
universes implies that the energy density of the created 
gravitons cannot increase with frequency if p -\- p > 0, 
i.e. if the DEC is not violated. Unfortunately, previous 
models which exploit this idea have a phase in their evo- 
lution where the matter entropy decreases. Specifically 
was considered a class of solutions which correspond to a 
viscous fluid with an equation of state given by p — —p. 
In this model the early phase (where the DEC is vio- 
lated) is smoothly connected to a radiation dominated 
evolution. Depending upon the sign of the bulk viscosity 
coefficient, the entropy of the sources driving the geom- 
etry can very well decrease 0, Q • 

Let us discuss the class of viscous cosmologies consid- 
ered in Ref. more in detail. In a flat FRW background 
the Einstein field equations in the presence of the bulk 
viscosity coefficient ^ can be written as 
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where the effective pressure is given by 
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In this case k = SttG, H — a/a, a{t) is the scale factor of 
the flat FRW metric and the overdot represents deriva- 
tion with respect to the cosmic time coordinate. In order 
to have the notation of the paper Q we must identify 
= 3/k and p' = P^ff. 
The equations Q|- © imply the energy balance 



3i7(p + Pe//)=0. 
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In order to have a cosmological model whose evolution 
violates the DEC only for a finite amount of time, in 
Ref. it is assumed that 



2H 
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This parametrization is very reasonable since the amount 
of violation of DEC is proportional to H. Effectively, 
from Eqs. and Q we have that for any solution 



P + Peff 
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and then a violation of the DEC implies that H > 0. 

In order to have a cosmology whose early phase (where 
the DEC is violated) is smoothly connected to a radia- 
tion dominated evolution, Giovannini considers the scale 
factor given by 



1/2 



(7) 



and then the self-consistent solution takes the form (see 
FigP) 
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One can immediately see by taking the limit t ±oo 
that C-oo(t) > 0, <e+oo(t) < and a±^{t) -> {±t)^^/\ 
So a{t) at the final phase of the whole evolution has the 
behavior of the radiation dominated era. 

Now, if we put the expression JSJ into Eq. @ we con- 
clude that the Giovannini parametrization implies that 
the local equilibrium pressure is given by p = —p [T(il |. So 
for t — > +00 the asymptotic solution to Eqs.Q and © 
is the following exact solution of the field equations 



a{t) = t 
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FIG. 1: We show the behavior of H, H and for viscous 
model O and We have cliosen ti = 1. We can see that 
H > (^ > 0) for t < (DEC is violated) and 7? < (^ < 0) 
for t > (DEC is preserved). 



Notice that Eqs. Q) and (O imply that 
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and then Peff p/3 for t +oo, while the local equi- 
librium pressure behaves always as p = —p. This is illus- 
trated in Fig|21 

Let us now consider some physical aspects of the dis- 
cussed viscous cosmological model. The thermodynam- 
ical entropy associated with the bulk viscosity © de- 
creases for t > 0. This implies that, for large positive 
cosmic time values, the second law of thermodynamics 
would be violated 0. However, as was discussed in 
Ref. 0, this statement might be not justified since it 
only takes into account the matter entropy but not the 
entropy of the geometry itself. This implies that, for a 
FRW background, in the framework of a well defined ex- 
tension of the second law of thermodynamics, one should 
include both the entropy connected with matter and the 
entropy connected with the FRW background. In this 
context, for positive cosmic time values, the decrease in 
the entropy of the sources may be compensated by the 
growth of the entropy of the FRW background. 

Unfortunately, at the present such an extension of the 
second law of thermodynamics is ambiguous even for cos- 
mological models which do not violate the DEC (see for 
example jTll|V So in this paper we are interested in find- 
ing viscous cosmological models which do not need such 
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FIG. 2: We show the behavior of the local equilibrium 
pressure p = —p (see Eq. @), efltective pressure Peff (see 
Eq. IIU 'I and the radiation pressure itp^ = l/U"^ for t > 0. 
We have chosen ti = 1. We can see that p is always negative 
and Peff behaves like p^^^ for t — » oo. 



an extension of the second law of thermodynamics. 

In the following we shall generalize the Giovannini class 
of solutions discussed above in order to include cosmo- 
logical scenarios which always have a positive viscous 
coefficient and thus satisfy the standard second law of 
thermodynamics. 

We will search for wider classes of solutions in the pres- 
ence of viscous sources with a linear barotropic equation 
of state p — for the local pressure. Specifically, we 
are interested in studying solutions which preserve the 
form of the scale factor (0 (in order to keep all physical 
properties explaining the slopes of the energy spectra of 
relic gravitons and the smooth transition to the radiation 
dominated era) but having a positive bulk viscous stress 
in order to have an increasing matter entropy during all 
evolution of the cosmic time. 

From Eqs. |^ and (0) we obtain that the bulk viscous 
coefficient for this kind of fluid may be written as 



2H 



(7 + l)H. 
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Notice that Eq. (|12|l implies that the parametrization (O 
has a state parameter 7 = — 1. So we can consider viscous 
cosmological models for which the state parameter 7 ^ 
— 1, in other words we shall find a self-consistent solution 
for the full set of Einstein field equations. The second 
term of (|12|l may be positive and then we can have a 
non-negative bulk viscosity for all cosmic evolution. We 
can see that this is possible for expanding universes (for 
which H > 0) with state parameter 7 > — 1. 

Now from the field equations JQl and (jSJl we have that 
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and then this new class of viscous models will have the 
same effective pressure (|ll|l . which behaves as Peff 
pjZ for t +00. 

The new class of solutions can be written as 



a(t) = h + + ^1 
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3(7 + l)v/t^+7f-4t 
^ = "^^^^ WT^) • ^ ^ 

From here we get that the viscous pressure is given by 



FIG. 3: Some curves for the bulk viscous stress given in 
Eq. are plotted. We set t\ = \ and 7 = —3/2, —1,0 as 
illustrative values. We can see that each curve has a phase 
where it is positive, and then it transits to a negative phase. 



_3^^. _ 4t-3(7+l)v/^ 

4(t2+t2)3/2 ' (1^) 



and then Peff takes the form of Eq. (|ll|l . 

Now from Eq. H15|l we can derive the general behavior 
of ^. Note that its numerator in general can be positive, 
negative, or zero. It can be shown that, if 
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FIG. 4: Some curves for K^{t) given in Eq. HSU are plotted. 
We set ti = 1 and 7 = —3, —7/3, 1/3, 1 as illustrative values. 



the bulk viscous stress is zero. This root is a real one 
if —7/3 < 7 < 1/3. This means that, in this range 
of the state parameter, the bulk viscous stress has a 
phase where ^ is positive (the matter entropy increases 
with time) and another phase where it is negative (the 
matter entropy decreases with time). In this case, if 
—7/3 < 7 < —1, the value of troot is negative and, if 
— 1 < 7 < 1/3, the value of troot is positive. This is shown 
in Fig|3| Note that the solution studied in Ref. 0,11] lies 
in this range so, for —7/3 < 7 < 1/3, we have a class of 
Giovannini-like solutions with 7 7^ —1. The solutions for 
which the matter entropy always increases with cosmic 
time lie out of the range —7/3 < 7 < 1/3. In this case 
the root H17|l does not exist and then the bulk viscosity 
is always negative or always positive. Effectively one can 
show that, for 7 < —7/3, the bulk viscous stress is al- 
ways negative, and for 7 > 1/3 the bulk viscous stress is 
always positive (see Fig. 

In conclusion, any viscous cosmology with a constant 
barotropic state parameter 7 > 1/3 will have a positive 
bulk viscous stress, and thus an increasing matter en- 
tropy during all cosmic evolution. This new class of solu- 
tions has the following asymptotic behavior at f — > +00: 
np 3/4^2, KFe// ^ 1/4*2, _^ (-3^ - l)/6t and al- 
ways p = 7p. From here we conclude that the radiative 
viscous fluid solution is more physically acceptable than 
other ones since, for it, at first order, we have Peff P, 
^ « and always p = p/Z. So, it should be emphasized 
that lit ^ +00 the scale factor expands as a — > t^/^ while 



the equation of state of the local equilibrium pressure p 
may be different from p/i. However, what counts is the 
effective pressure for which we have that Peff —>■ p/3 as 
it should. 

Notice also that the sign of the effective enthalpy, i.e. 
p + Peff, plays an important role. Effectively, we can 
see that the DEC is associated with Eq. ©, which im- 
plies that for iJ > DEC is violated, and for < 
DEC is preserved. On the other hand, as we have stated 
above, the slopes of the energy spectra of relic gravi- 
tons are crucially related to the sign of the effective en- 
thalpy. The energy density of the created gravitons in- 
creases with frequency ii p + Peff < 0, i.e. if the DEC is 
violated {H > 0). So the considered here viscous cosmo- 
logical models with a constant barotropic state param- 
eter 7 > 1/3, violate the DEC only for a finite amount 
of time. In this early stage {t < 0) the relic gravitons 
are created, and when the DEC is restored {t > 0) the 
universe exits to the standard radiation dominated stage. 
During all evolution of the cosmic time the bulk viscosity 
is positive. In this way we have shown that the possible 
bulk viscous influence on the slopes of energy spectra of 
relic gravitons makes sense even without violation of the 
standard second law of thermodynamics. 

Lastly, notice that we have considered the same form 
for the background metric 0; thus all calculations re- 
ported in Ref. ^ are preserved since they rely mostly on 
the specific time dependence of the scale factor, rather 
than on the bulk viscosity coefficient. Thus we come 
to the same conclusions concerning the amplification in- 
duced by the background 10 in the proper amplitude of 
the gravitational waves reported in Ref. |^. It is im- 
portant to stress here that these results imply that the 
existence of tilted spectra of relic gravitons can be con- 
nected, in the framework of general relativity, with the 
violation of the DEC induced by a bulk viscosity com- 
patible with the standard second law of thermodynam- 
ics. Finally, going back to the original Giovannini mo- 
tivation 0, we conclude that these results indicate the 
possibility of the existence of growing spectra of gravita- 
tional waves at frequencies of the interferometric devices 
which can put bounds on the possible violation of the 
DEC occurring in the early universe. It is interesting to 
note that growing energy spectra can also be obtained in 
the framework of non-Einsteinian theories [3, [131 . 
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